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ABSTRACT 

Context. Magnetic twist is thought to play an important role in coronal loops. The effects of magnetic twist on stable magnetohydro- 

dynamic (MHD) waves is poorly understood because they are seldom studied for relevant cases. 

Aims. The goal of this work is to study the fingerprints of magnetic twist on stable transverse kink oscillations. 

Methods. We numerically calculated the eigenmodes of propagating and standing MHD waves for a model of a loop with magnetic 

twist. The azimuthal component of the magnetic field was assumed to be small in comparison to the longitudinal component. We did 

not consider resonantly damped modes or kink instabilities in our analysis. 

Results. For a nonconstant twist the frequencies of the MHD wave modes are split, which has important consequences for standing 
waves. This is different from the degenerated situation for equilibrium models with constant twist, which are characterised by an 
azimuthal component of the magnetic field that linearly increases with the radial coordinate. 

Conclusions. In the presence of twist standing kink solutions are characterised by a change in polarisation of the transverse displace- 
ment along the tube. For weak twist, and in the thin tube approximation, the frequency of standing modes is unaltered and the tube 
oscillates at the kink speed of the corresponding straight tube. The change in polarisation is linearly proportional to the degree of 
twist. This has implications with regard to observations of kink modes, since the detection of this variation in polarisation can be used 
as an indirect method to estimate the twist in oscillating loops. 
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1. Introduction 

It is anticipated that magnetic twist plays a considerable role 
in the structure of coronal loops. The highly dynamic photo- 
sphere and chromosphere can introduce twist in the magnetic 
field of the solar corona. For example, footpoint rotation and 
shear motion may result in twisting of coronal loops (see e.g., 
iBrown et al.l l2003b . A newly emerged magnetic field is also 
supposed to be twisted during the buoyant evolution through 
the convection zone, meaning that loops may emerge already 
twisted (see e.g., Moreno-Insertis & Emonet 1996; Ho od et al.l 
|2009|) . Twist also considerably contributes to the eruption of 
many prominences and CME s (coronal mass ejections) (see e.g., 
iHood & Priestlll979HPriest e t al. 1989). 

In the theoretical analysis of magnetohydrodynamic (MHD) 
waves, the effects of a curved loop axis are usually omitted. The 
stability of straight tubes with t wisted magnetic field has be en in- 
vestigated by e .g., IShafranovl (Il957h . iKruskal et all ( Il958l) . and 
ISuydamI (Il958h . The role of line tying c onditions in th e context 
of coronal loops has been studied by e.g.. lRaadul (Il972h while the 
effe ct of gas pr essure on the kink instability has been analysed by 
e.g.. lGiachetti et al. (1977). The combinati on of both eff e cts h as 
been stud ied by Hood & Priest (1979), Hood & Priest' (T981), 
iHood et al. (1982),.Einaudi & va n Hoven tl983.) . and Velli et al. i 
(1990). More recently. Zaqarash vifi et al.l ( 120101) and lDiaz et al.l 
(201 ll) have considered the stability of twisted magnetic flux 
tubes with mass flows along the field lines. The main aim of 
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the previous works is determining the instability threshold and 
calculating the growth rate of kink unstable modes. 

Stable propagating MHD waves have been studied previ- 
ously, but almost always under particular conditions that some- 
how reduce the mathematical complexity. Often the MHD waves 
are incompressible and/or the azimuthal component of the mag- 
netic field varies linearly with distance up to a certain po- 
sition where it drops to zero or all the way up to infinity. 
The first choice leads to a surface current while the second 
choice causes magnetic energy to d iverg e. Early studies were 
performed by iGoossens et alJ d 19921) and iBennett et all (119991) 
in the incompressible regime. The analysis of sausage modes 
(m = 0) has been carried out bv lErdelvi & FedunI (I2006ll2007h . 
Kink modes assuming compressibl e motions have be e n con - 
sidered by Erd elvi & Cartel (l200d) : ICarter & Erdiivil (l2008l) : 
lErdelvi & Pedum (|201(I) . In these works standing transverse 
waves have not been considered and only propagating waves 
have been analysed. 

In particular, 'RudermanJ (l2007h (see also review by 
Ruderman & Erdelyi 2009) studied standing waves and con- 
cluded that twist does not af fect kink mode s, only fluting modes 
are modified (see also Gooss ens et al.ll992h . This result is due to 
the particular choice of constant magnetic twist, i.e., twist is in- 
dependent of the radial coordinate. The conclusions of the work 
of Ruderman (2007) have some authors to regard twist as unim- 
portant for transverse kink oscillations. 

In this paper we study the effect of an azimuthal magnetic 
field component on the frequencies of MHD waves for an equi- 
librium model with nonconstant magnetic twist. We are espe- 
cially interested in the effect of twist on standing kink modes, a 
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problem that, to our knowledge, has not been addressed in the 
context of transverse coronal oscillations (with the exception 
of standing fluting modes, studied in Ruderman 2007). Thus, 
our focus is on the fingerprints of twist on transverse kink os- 
cillations but in the stable regime. This might be relevant for 
interpreting kink oscillations, from the early transvers e kink 
oscillations reported with TRACE (lAschwanden et all Il999t 
[Nakariakov et al. 1999 ), to the last AIA observations of trans- 
verse oscillating loops dAschwanden & Schriivedl201 lb . 

2. Equilibrium model of the magnetically twisted 
loop 

The model for the coronal loop that we use in the present paper 
is a straight cylinder with constant density inside the cylinder 
and constant density outside the cylinder. Thus, the coronal loop 
is modelled as a density enhancement characterised by a den- 
sity contrast pi/pe with pi being the internal density and pe the 
external density. 



Pi, 0<r<R, 
Pe, r > R. 



(1) 



Atr-R the density changes in a discontinuous manner and R is 
the radius of the loop. 

Under coronal conditions it is a good approximation to ne- 
glect plasma pressure compared to magnetic pressure. This zero- 
j6 approximation removes the slow magnetoacoustic waves from 
the analysis. In the present paper we allow the magnetic twist to 
be nonconstant. The model starts from the equation of magneto- 
static equilibrium 



dr r 
where 



= bI + bI- 



(2) 



(3) 



Here the azimuthal component of the magnetic field is pre- 
scribed and Eq. ^ is solved for B-(x) with x - r/R. The solution 
is 



B^Jx) = B^M - Bl(x) 



Jo 



-ds. 



The azimuthal component of the field is written as 

B^(x) = A fix). 



(4) 



(5) 



where A has the dimension of a magnetic field and f(x) is a di- 
mensionless function satisfying the condition /(O) - or better 
f(x) — > for X — > 0. If we insert Eq. (|5]l in Eq. (01 we find 



bI(x) = B^(0) - A^G(x), 



G(x) = fix) + 2 



Jo s 



ds. 



(6) 
(7) 



Notice that G(0) = and G{x) > 0. 

A particular choice for f{x) used in the present paper is 

f 0, Q<x<p, 
f{x) ^ < {x- p){q- x), p <x<q, 
10, x> q. 



(8) 



where p and q are constants that we can freely choose. The 
choice /5 = is an option and that ^ > 1 is also allowed. For 



this model we have the freedom to have twist inside and also 
outside the tube. The function fix) attains its maximal value 
fmax = iq- p)^/^ at x„, ^(p + q)/2. Hence 



B^ix) = aBM-^. 

Jmax 



(9) 



In order to relate the strength of the azimuthal component B^ 
with respect to the dominant longitudinal component, we intro- 
duce a: 



B,r, 



B^iO) 



(10) 



hence 

A = -^BM- 

Jmax 



(11) 



For the function f{x) specified by Eq. ^ it follows that 

f 0, 0<x< p, 
G{x) = I > 0, p < x<q, 
[ G(q), x> q. 

This function is calculated using Eq. O. 

In summary, the equilibrium magnetic field is defined as 

B, — Bq, B^p — 0, 

for X < p, 



(12) 



(13) 



B, = ^Bl - A^G(x), B^ = Aix - p)(x - q), (14) 
for p < X < q, with 



G{x) -(x- p)^(x - q)^ + c^{x^ - p^) + c-i{x^ - p^) 
+C2ix^ - p^) + ciix- p) + ci (ln(x) - \n(p)) , 



ci = -4pq{p + q), 

C2 ^ (p + qf + 2pq, 
4 

C3 = --^(p + q). 



1 



C4 



o 2 2 

ci ^ 2p q , 
and 



(15) 

(16) 
(17) 

(18) 

(19) 
(20) 



B^ = 0, B, = ^Bl-A^G(q), (21) 
for X > q. 

3. MHD Eigenmodes 

3.1. MHD equations and numerical method 

The motions superimposed on the equilibrium model given in 
the previous section are governed by the linearised ideal MHD 
equations in cylindrical coordinates. These equations can be 
found in the appendix. Several terms involving radial derivatives 
of the equilibrium magnetic field are present in the equations. 

In the MHD equations we have performed a Fourier analysis 
in time, with <jj the frequency, and also in the azimuthal direc- 
tion, m being the azimuthal wavenumber. For the longitudinal 
direction (along the tube axis) two situations are studied. The 
first situation corresponds to fully propagating waves, allowing 
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Fourier decomposition in z, k being the longitudinal wavenum- 
ber. In this case the equations are solved in the radial coordi- 
nate r. The second situation corresponds to standing waves with 
line-tying conditions at the footpoints of the loop. This prob- 
lem does not allow Fourier analysis in z, and the equations are 
solved in 2D (in r and z). For the analysis of propagating waves 
the system of equations reduces in the ideal regime to the Hain- 
Liist equation. In general, it is not possible to find analytical so- 
lutions to this equation when the twist profile is nonconstant. 
For this reason the MHD equations are numerically solved using 
the code PDE2D (Sewell 2005) in the two situations (ID and 
2D) explained above. More det ails about the numerical method 
can be found in Terra daFet al.l ([2006). In the present paper we 
avoided possible Alfvenic resonances on purpose, interesting as 
they might be. In addition, we stayed away from a strong az- 
imuthal field that can cause kink instabilities of the magnetic 
flux tube. This means that the eigenfrequency w is a purely real 
magnitude in our analysis. 

3.2. Effect of the azimuthal magnetic field on MHD waves 

Since we adopted the zero-/? approximation, the only character- 
istic frequency in the system is the local Alfven frequency cja- 
The effect that the azimuthal component of the magnetic field 
has on the MHD waves can be appreciated from the effect of B^ 
on 0?.. The square of the Alfven frequency is defined as 



1 /m \2 
i-B^ + kBA 

Hop \ r I 



This expression can be written as 

Bl 

The function <t>(r) is 



3? 1 ,/ m \2 
^-^(kRf(l + -^(r)] 
()PR^ \ n I 



<D(r) 



kR{rlR)B- rB, 



(22) 



(23) 



(24) 



where we have defined L that k - njL. <b{r) is the twist of the 
magnetic field lines and it is 27r times the number of windings 
of the field around the loop axis over a distance L. Equation (l23T l 
suggests that even a relatively small y?-component of the mag- 
netic field in the sense that \B^\ « \B,\, might have a relatively 
important effect on the MHD waves. The relative importance of 
B^p to B~ is controlled by the amount of twist of the magnetic 
field, with the length L being a relevant factor 

The ratio of the contribution to the Alfven frequency due 
to the azimuthal component to that due to the longitudinal 
component is 



B. 



kR (r/R)B, 



(25) 



We computed the maximal value of the ratio for the model dis- 
cussed in Section |2] We considered a weak azimuthal field and 
a relatively low value of << 1. To first order in a (small a 
approximation) 



B. ^ 5.(0), 

B^ = aBMm/fm 

Hence 



B^ 



ir/R)B, 



m 1 

J ma. 



(26) 
(27) 



(28) 



For the prescription of f{x) given by Eq. ([8]) we already know the 
value of x,„ and fmax- In addition, the function /(^)/x is maximal 
at jc» - -^/pq, and f(xt)/xt = p + q-2 s/pq. Hence the maximal 
value of 



Bu) 



{rlR)B, 



- a- 



(q - py 



{p + q-2 Vm) ■ 



(29) 



For example, for the values p - Ijl and q - ill used in our 
calculations we find that 



Bu) 



kR {rlR)B, 



1.08 m—. 
kR 



(30) 



As an example we consider a loop that is 10 times longer than 
the radius. For the fundamental mode along the tube axis it fol- 
lows that kR - nllQ. So for a - 1/10 the ratio is about 0.34. 
Thus, although the azimuthal component of the magnetic field is 
weak, its effect on the Alfven frequency is substantial. For real- 
istic loops it is a competition between the two quantities kR and 
a which determine the ultimate effect. 

We omitted resonantly damped eigenmodes. This choice re- 
stricts the maximal twist allowed in the system or conversely the 
allowed combinations of a and kR. Even for a constant density 
profile like the one considered in this work (see Eq. ([T]i), wa is 
a function of position because B^ varies with position and the 
corresponding term is multiplied by the factor 1 jr. A necessary 
(but not sufficient) condition for the absence of resonances is that 
maxw?. < min . This conditions is equivalent to 



max <b <ni^ y/pi/Pe ~ 1 ) ■ 



(31) 



For pi/pe = 3 this implies that max (£> < 0.737r. 

Condition (13X1 1 can be rewritten in terms of the function f(x) 

as 



^ fmax 
< X» 

kR fix,) 



( VPi/Pe - l)- 



(32) 



For the typical values that we use in the computations the in- 
equality is a/kR < 0.67. So we can define the critical value of a 

as 



aj = kRx,jj^i^yjpilps - l). 



fM 



(33) 



Then for values a > aj the azimuthal component of the mag- 
netic field introduces resonances even in configurations with 
constant density. In addition to ffy we can consider the situation 
when the azimuthal component dominates the variation of 
<x>\. For m - 1 this happens when max(I>(r) > n (see Eq. (l25ll). 
For the function f{x) defined in Eq. dS) this can be reformulated 
into 



kR 



> X, 



fn 



fM 



This forces us to define a second critical value for a 

fmax 



Qm — kR Xt 



f(x,y 



(34) 



(35) 



which differers from ffy by the factor ( yJpi/Pe - l)- For the 
present model with p = 1/2 and q = 3/2, - l.09kR. Below 
we used sufficiently low values of a compared to ffy and um. 
This implies that even in the thin tube approximation we can use 
the small a approximation. 
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4. Results: propagating waves 

We start by solving the eigenvalue problem for propagating 
waves using different pairs of the wavenumbers m and k, which 
are assumed to be positive. We use a particular choice of param- 
eters, pi/pe - ^, p - ^ /2, q - 3/2. In Figure[T]the dependence 
of the eigenfrequency on a, which measures the amount of twist, 
is plotted. For the pair im,ki) the frequency is w = kiC]^ when 
a - 0, i.e., we recover the known result in the absence of twist 
(in the thin tube approximation), which is the kink speed 



Ck = 



\ 1 +Pe/Pi 



VAi- 



(36) 



where VAi = Bo/ VjUoPi is the internal Alfven speed for a straight 
field. 

Figure [1] shows that the frequency increases monotonically 
with a (see continuous thick line). Exactly the same curve is 
obtained for the wavenumbers (-m,-ki), meaning that the re- 
sults are degenerate with respect to a change in sign of the two 
wavenumbers. In constrast, for {-m,ki) and (m,-k\) the fre- 
quency decreases with a (see dashed thick line). This behaviour 
reflects the fact that, in contrast to the purely vertical magnetic 
field model, twist breaks the symmetry with respect to the propa- 
gation direction when the sign of one of the two wavenumbers is 
changed (see also the results of Vasheghani Farahan i et al.ll2010i 
for torsional Alfven waves). The increase in frequency of the 
modes (m, k\) and (-m, -ki) with respect to the un-twisted case 
is around 4% while the decrease for {-m, k\) and (m, -k\) is 8% 
for a = 0.01. This means that a very weak twist (of about 1%) 
produces an effect on the frequency that is not negligible. The 
maximum value of a in our calculations has been chosen to 
satisfy that a < aj and a < au for the reasons explained in 
Section[l2] 

Figure [1] also suggests that a linear variation of w with a is a 
good approximation of the actual dependency. 



ii) — kiC)^ + aa, for (m,ki), {—m,—ki), 
to ^ kicy_ — aa, for (m,—ki), (—m,ki). 



(37) 
(38) 



a being the sl ope of the linear approx imat ion. The tw ist profiles 
considered bv lGoossens et alJ (Il992h and lRudermanI (l2007l) are 
special since for those cases a = 0, i.e., there is no frequency 
change for different wavenumbers (for example, m = ±1). 
However, for nonconstant twist configurations, like those studied 
here, there is always a frequency split. 

Similar behaviour occurs for different longitudinal wave 
numbers, the curves are just shifted vertically (see continuous 
and dashed thin lines). The slope of the curves is independent 
of the value of k, as we can see in Figure [1] (compare the curves 
for ki and k2), but as we show below, it depends on the par- 
ticular twist profile. According to Figure[T]for a fixed frequency 
and fixed value of twist, there are always four possible solutions, 
(m, ki), (-m, -ki), (-m, k2), (m, -k2) (see the intersection of the 
solid thick and dashed thin lines around oj = 7.85 x 10 ' and 
a = 3.5 X lO^^*). This property will be used as the basis to con- 
struct approximate standing solutions in the following sections. 

A qualitatively similar behaviour occurs when the parame- 
ters are changed. In Figure |2] the results are plotted for p - Q 
and q = 2 (see thin lines). Now the slope of the curves is steeper 
than for the case p - Ijl and ^ = 3/2 but the curves still show 
a linear behaviour with a. Since for the case p - and q - 2 
the region where the twist is localised is larger, the effect on the 
frequency is more pronounced. Indeed, we tested different non- 
constant twist profiles and the curves are always linear with a or 
equivalently with (for << |Bj:l). 



0.084 



0.082 



0.080 



■S' 0.078 
3 



0.076 



0.074 



0.072 




(—m, k^) ( ra, -fc^) 



0.070 I I I I I I I I I 

0.000 0.002 0.004 0.006 0.008 0.010 



Fig. 1. Frequency as a function of twist {a). The curves are obtained 
by solving the eigenvalue problem in ID. The dotted line represents the 
solution for the untwisted case. For this particular example, pjp^ = 3, 
p = 1/2, q = 3/2. The two sets of solutions are associated to m = 1 and 
ki = n/50R (thick lines), and k2 = 7t/4SR (thin lines). 



We now aim to know the effect of twist on the eigenfunc- 
tions. It can be shown from the MHD equations (see Appendix) 
that if Fourier analysis is performed in both the azimuthal and 
longitudinal direction, like in the present case, then there is a 
phase shift between the radial component (pure real function) 
and the azimuthal and longitudinal components (purely imag- 
inary functions). In Figure [3] the radial dependence of vv, 
and Vj is plotted for the modes (m,k) (continuous line) and 
(-m, k) (dashed line). The same figure is obtained for the modes 
(-m, -k) and (m, -k). We see that the radial dependence of the 
modes is essentially the same inside the tube and only in the ex- 
ternal medium near the tube boundary the dependence is slightly 
different (compare also with the dotted line, which represents the 
eigenfunction for the classical untwisted tube). The azimuthal 
velocity component shows more relevant differences for the two 
modes. The mode (m, k) has a reduced shear at the tube boundary 
while the opposite behaviour is displayed by the mode {-m, k), 
showing an increased shear with respect to the un-twisted case. 
This behaviour is also present in the v- component, which is in- 
troduced by the presence of twist. Since we are in the zero-/? 
approximation, the velocity along the magnetic field is zero, i.e., 
V|| = V ■ Bo/|Bol = 0, as expected. The v, component has a jump 
at the tube boundary and is different from zero only in the region 
where there is twist (in this particular example between r = 1 /2R 
and r = 3/2R). Notice that the amplitude of this velocity com- 
ponent is rather small in comparison with and since we are 
considering a weak twist. 
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Fig. 2. Frequency as a function of twist {a). For this particular exam- 
ple, Pi/pc = 3, p = 1/2, q = 3/2 represented with thick lines while 
the case p = 0, q = 2 is plotted with thin lines. The dotted line repre- 
sents the solution for the untwisted case. The horizontal continuous line 
is obtained by numerically solving the eigenvalue problem in 2D. The 
same notation as in Figure[T]is used. 

5. Results: standing waves 

So far we focused on both azimuthally and longitudinally prop- 
agating waves. Our interest is now on the analysis of standing 
waves, if they exist in the twisted configuration. 

5.1. Standing in (p, propagating in z 

An azimuthally propagating wave, for example the mode m- 1, 
produces a motion that displaces the whole tube and its axis is 
moving following a circular path because of the propagating na- 
ture of the wave in ip. If the motion is clockwise for m then it 
is anti-clockwise for -m. For the tube to oscillate transversally, 
the propagating m mode has to be combined with the -m mode. 
When there is no twist, the frequency of the mode (m, k) is the 
same as the for the mode {-m,k) (oj = kc^ in the thin tube 
limit), and the superposition is trivial. In the presence of twist 
the situation is more complicated since, as we have shown in 
Section |4] the modes (m, k) and {-m, k) have different frequen- 
cies. This means that with a single longitudinal wavenumber it 
is not possible to have a standing solution in the azimuthal direc- 
tion. However, we have shown that we can always find another 
pair of wavenumbers with the same frequency. This means that 
the linear combination of the modes {m,ki) and (-m,k2) will 
lead to the standing solution in ip we are looking for. Formally, 
we have to combine the two waves, 

y - fif) 6xp / {(jL>t + imp + kiz) + g{r) exp / {ut - imp + k2z), (39) 

where y is, for example, any of the velocity components of 
the eigenfunction. Note that f{r) and g(r) are different func- 
tions since they are associated to different pairs of wavenumbers 




2.0^ 
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Fig. 3. Dependence of the three velocity components of the eigenfunc- 
tion with the radial coordinate. For this particular example, pjp^ - 3, 
p = 1/2, q = 3/2 and a = 0.01. The continuous and dashed lines corre- 
spond to the modes (m, k) and {-m, k), while the dotted line represents 
the mode (m, k) in the absence of twist (a = 0). 



{m,ki) and (-m,k2). However, in the thin tube approximation, 
when k\R <s: 1 and k2R <K 1, the two functions are essentially 
the same (/(r) a: g{r) or/(r) a: -g{r) as Figure[3]suggests). This 
is a nontrivial assumption, which might not be a good approxi- 
mation when we relax the thin tube and the small a approxima- 
tions. Using this approximation, Eq. ( |39] | reduces to 



y - fir) COS {nup + kjz) exp / (cot + kz). 



(40) 



where we have introduced for convenience the following nota- 
tion 



k = 
kj — 



k\ + k2 
2 ' 
ki - k2 



(41) 
(42) 



If instead of the sum of the modes (m, ki) and (-m, k2) we choose 
the difference as the linear combination, we find 



y - .fi'') sin {m(p + kjz) exp / ((Lit + kz). 



(43) 



It is straightforward to construct the standing solution in ip but 
propagating in the positive z-direction (note that Eq. (l40l i repre- 
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sents a wave propagating in the negative z-direction). The solu- 
tion is based on the superposition of the modes {-m, -ki) and 
(m, -k2), and it is 



y = /('■) cos (mip + kjz) exp / (cut - kz). 



(44) 



Again the sinus solution is obtained when we subtract the modes. 

Equations ( l40l l and ( l44l l represent propagating waves in the 
z-direction and produce a transverse displacement of the tube so 
that now it depends on z through the term kjz. Thus, the direc- 
tion of the transverse oscillation of the loop varies as the wave 
propagates along the tube. This represents a change in the polar- 
isation of the transverse motion and is one of the main effects of 
twist on the eigenmodes. For the untwisted tube the polarisation 
of oscillation is always the same since kj - 0. 

The choice of and kz is not arbitrary since the modes with 
the pairs of wavenumbers {m,ki) and (-m, ^2) (the same ap- 
plies to the combination {-m,-k\) and (m, -^2)) need to have 
the same frequency. The appropriate wavenumbers can always 
be determined numerically for a given frequency, as Figure [1] 
suggests. However, if a is small, this selection is straightforward 
because we can use the approximations for the frequency (see 
Eqs. (137])-(138])) to calculate ky and ki. 



o) + aa 



Ck 



(45) 
(46) 



where a is the slope of the curves. Now using Eqs. ( l4TT i and ( l42l i 
we have 



k 
kf 



a a 
Ck ' 



(47) 
(48) 



From the first equation we conclude that the dispersion relation 
of the wave in the twisted case, u - kc^, is exactly the same 
as in the untwisted tube (under the thin-tube approximation and 
for weak twist). Thus, for weak twist, there is no eff'ect on the 
frequency of the propagating waves in z and standing waves in 
(p. Nevertheless, from the second equation we find that the po- 
larisation of the motion shows a dependence with twist (kj + 0). 
This dependence is linear with a. For the un-twisted case (a = 0) 
we recover the well-known solution of the transverse kink mode 
without variation of polarisation along the tube. 

5.2. Standing in z, propagating in ip 

Now we seek solutions that are standing in the z-direction. 
These solutions have to satisfy line-tying conditions, i.e., the 
three velocity components must be zero at the footpoints of the 
loop, located at z = and z - L, where L is the total length of 
the tube. We proceed as in the previous section and select the 
appropriate superposition of modes. A suitable choice is {m,k\) 
and (m, -kz), 

y - f(f) (exp i {cot + imp + kiz) + exp / {cjt + nup - kzz)) , (49) 
which can be written as 



y - f(r) cos (kz) exp ; (ojt + imp + kjz). 



(50) 



The other choice corresponds to the modes (-m, -ki) and 
(-m, +k2), which represent a propagating wave in the positive 
^-direction. 



As in the previous expressions, a sinus is obtained instead of a 
cosinus when the difference of the modes is considered instead 
of the sum. 



y - fir) sin (kz) exp / (cjt + nup + kjz), 
and 

y - /('■) sin (kz) exp / (cjt - imp - kjz). 



(52) 



(53) 



Strictly speaking these modes cannot be classified as fully stand- 
ing modes, because of the z dependence in the exponential. 

Applying line-tying conditions imposes restrictions on the 
wavenumbers. Using the sinus solution given by Eq. ( l52l i or (|53] ) 
allows us to easily impose that the velocity must be zero z = 
and z - L. The first condition is trivially satisfied while the 
second leads to 



sin (kL) - 0, 
and thus 



k — —n, n — 1, 2, 3 . 
L 



(54) 



(55) 



y = fi^f) cos (kz) exp / (cot - nup - kjz). 



(51) 



This is essentially the same situation as in the untwisted tube, 
we have a discrete set of wavenumbers satisfying the boundary 
conditions. 

In constructing the standing solution in z we have assumed 
that this solution is the superposition of two propagating waves 
with the same frequency, different longitudinal wavenumbers 
and the same radial dependence. As we have explained, this is 
only true within limit kR <g: 1 . In principle, we cannot assume 
a Fourier analysis in the z-direction because of line-tying con- 
ditions and consider single modes if we aim to obtain the full 
solution to the problem. For this reason, we now calculate the 
eigenmodes by solving the problem in r and z for a fixed in. This 
will allow us to compare the results with the semi-analytical re- 
sults given by Eqs. (l47b and ( |48] ). We used the two-dimensional 
version of the code PDE2D to calculate the modes with line- 
tying conditions at z = and z - L. In Figure |2] the eigenfre- 
quency is plotted as a function of twist. Clearly, the frequency 
is almost independent of a (see the continuous horizontal line), 
which agrees with the analytical results obtained within the limit 
of small a given by Eq. (l47l l (compare also the continuous hor- 
izontal line with the dotted line, which corresponds to the un- 
twisted case). Therefore, it is a good approximation (within the 
limit of small a) to assume that we can construct standing eigen- 
modes by summing two propagating waves with the same fre- 
quency but slightly diff'erent wavenumbers. 

A plot of the eigenfunction at r = and 1^ = as a function 
of z is shown in Figure |4] It is worth to mention that for the 2D 
problem the eigenfunction is always a complex number (this can 
be seen from Eqs. (IA.ll i- (IA.6b in the appendix) although the fre- 
quency is still a real magnitude since we do not consider unstable 
modes or resonances. For this reason we represent the modulus 
and the phase (6) of the eigenfunction in Figure |4] The modulus 
has a sinusoidal profile while the phase shows a linear depen- 
dence with z. This is exactly the behaviour predicted by Eq. (l52T i. 
The slope of the phase is precisely the magnitude kj, which is 
negative according to Eq. ( |48] ). The in = -1 mode shows the 
expected positive slope, in agreement with Eqs. (|53] | and (l48b . 

5.3. Standing in ip and z 

Once we have derived the solutions that are standing, or equiv- 
alent to standing, in one of the directions, it is straightforward 
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inside the tube f(r) ^ g{r) - const in the thin-tube approxima- 
tion, and the displacement in z has a node at the axis. Thus, in 
Cartesian coordinates at the axis {x - y — 0) 



= ^0 cos (^tZ + Vo) sin (fe) sin (wf) > 
= ^0 sin {kjz + (fio) sin (kz) sin (cot) , 



(60) 
(61) 
(62) 



^0 being an arbitrary ampHtude and ipo an arbitrary phase. These 
expressions show that the tube does not move along the axis 
while the displacement in the plane perpendicular to the axis 
is given by the product of two sinusoidal terms, one involving a 
scale associated to the twist, kj, and another related to the longi- 
tudinal wavenumber, k. At the axis 



0.10 




-0 10 \ I I I I 1 

10 20 30 40 50 



z/R 

Fig. 4. Eigenfunction as a function of z at the loop axis calculated nu- 
merically using line-tying conditions. The modulus and phase are plot- 
ted for the modes m = ±1. 

to address the full standing problem. There are two ways to pro- 
ceed that must be equivalent, either we add the solutions that 
are standing in (/? but propagating in z, or we combine the solu- 
tions that are (almost) standing in z but propagating in if. For ex- 
ample, from the combination of the solution given by Eq. ( l40l l. 
representing a standing wave in if but propagating in z, minus 
the complementary wave travelling in the opposite z-direction 
given by Eq. (l44l i. we find the "full" standing solution 

y - fi'') cos (mtp + kjz) sin (kz) exp (iait). (56) 

It is easy to check that we obtain exactly the same expression 
by combining the solutions given in Eqs. (l52]i-([53Tl for mainly 
representing standing waves in z but propagating in if. 

Based on Eq. (l56l l and using the proper phase shifts between 
the components of the displacement, we obtain the following 
eigenfunction 

& - /('■) cos (mif + kjz) sin (kz) sin (cot) , (57) 

- sir) sin (mip + kjz) sin (kz) sin (cot) , (58) 

- h(r) sin (mtf + kjz) sin (kz) sin (cut) . (59) 

Similar expressions have been used in the past as approximate 
solutions to test the stability of different twisted magnetic con- 
figurations. However, as far as we know, these expressions have 
not been applied to stable kink waves. 

From Eqs. (I57b-(l59b we can easily calculate the displacement 
of the tube axis, which is a relevant quantity when it comes to 
observations of transverse loop waves. According to Figure [3l 



+ ^2 + ^2 ^ ^0 sin (kz) sin (cat) , (63) 

meaning that the total displacement is exactly the same as in 
the un-twisted case. The difference is in the components of the 
displacement which have a different weight depending on the 
position along the tube. 

In Figure |5] the displacement of the tube axis is plotted for 
a situation with moderate twist. Although the calculations in the 
previous sections were made very weak twist, we expect that 
for higher values of the component the motion of the axis is 
qualitatively the same. The differences in the polarisation of the 
motion with respect to the untwisted situation (see dotted line) 
are evident in this plot. The motion of the tube is even more clear 
in Figure |6] where the displacement of the tube is represented in 
three-dimensions at two different times during the oscillation. 
We clearly see the departures from the purely sinusoidal dis- 
placement for the untwisted case. In the top panel of Figure |6] 
the displacement is closer to the first longitudinal harmonic than 
to the fundamental mode. However, at later stages of the oscil- 
lation, see bottom panel of Figure |6] the motion resembles the 
fundamental mode. Thus, twist can have a pronounced effect on 
the displacement of the tube axis during the oscillation. 



1.0 




-0.5 L , I ^ , . , . ^ ^ I 

0.0 0.2 0.4 0.6 0.8 1.0 

z/L 

Fig. 5. Displacement of the tube axis as a function of position along 
the twisted tube. The continuous line corresponds to the displacement 
in the x-direction while the dashed line represents the displacement in 
the y-direction. The dotted line is the purely sinusoidal displacement in 
the X- direction for the untwisted case. In this example, = -0.75k, 
ifo = 0, where ^o/^ = 1- 
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Fig. 6. Snapshots of the displacement of the tube with twisted mag- 
netic field at two different times during one oscillation period. In this 
example, kj = -0.75fc. The x-axis is horizonal, the y-axis is vertical 
and the z-axis is perpendicular to the xj-plane. Line-tying conditions 
are applied at the front and back planes (z = and z = L). 



6. Conclusions and discussion 

From inspection of the local Alfven frequency it might be an- 
ticipated that an azimuthal magnetic field can have a relevant 
effect on MHD waves even if \B^\ « It is the competi- 
tion between the two quantities kR and a - B^ i„ax/B,(Q) that 
determines the extent of the effect. In this work the quantity a 
is a small quantity by choice since we aim to focus on magnetic 
fields that are dominated by the longitudinal components. The 
quantity kR is often small as most of the loops are much longer 
than wide or the longitudinal wavelength is much longer that the 
radius of the loop. So it is really a competition between two small 
quantities kR and a. In the present investigation we have delib- 
erately excluded the regime a/kR ~ 1 and focused our attention 
on a < a J and a < au- In this way we avoid the interesting 
complications of resonant absorption. However, restricting our- 



selves to these values of a means that the azimuthal magnetic 
fields are in a sense very weak, and that the effects on the MHD 
waves by azimuthal magnetic field are less pronounced than they 
effectively are. 

In general, the frequencies of MHD waves in models with 
nonconstant magnetic twist are nondegenerate compared to their 
counterparts in models with constant magnetic twist. This means 
that the frequencies of propagating waves with wavenumbers 
(m, K) and (-m, -K) are equal but different from the frequency 
of the modes (-m, k) and (m, -k). Twist breaks the symme- 
try with respect to the propagation direction. This property has 
consequences for standing modes that require the superposi- 
tion of propagating waves. We have shown, using an approxi- 
mate method valid for weak twist, that the frequency of standing 
waves is unaffected by the presence of twist. Under such con- 
ditions the tube oscillates transversally at the characteristic kink 
frequency. 

However, the most interesting effect of magnetic twist on 
standing oscillations is in the change of polarisation of the trans- 
verse motion along the tube. We have analytically demonstrated 
that the change in the direction of the transverse displacement 
of the tube is linearly proportional to the twist. Moreover, the 
polarisation of the displacement of the axis along the tube is 
simply described by the product of two sinusoidal terms in- 
volving a wavenumber associated to the twist, and a longitu- 
dinal wavenumber that satisfies the line-tying boundary condi- 
tions at the footpoints of the loop (this is very similar to the sit- 
uation for a tube wi th a siphon flow along the axis studied by 
iTerradas et al.ll201lb . For a very weak twist, for instance 1/10 
windings of the field along the loop axis, twist apparently intro- 
duces a motion perpendicular to the dominant direction of oscil- 
lation that is around nine times smaller in amplitude. The effect 
on the polarisation of the motion is small because the twist in the 
tube is very weak in this particular example. When the twist in- 
creases, the change in polarisation becomes stronger, as is shown 
in Figures |5] and |6] An important result here is that a weak twist 
can produce displacements in any direction perpendicular to the 
unperturbed tube axis. Thus, in contrast to the effect of curvature 
of the magnetic field lines and the effect of nonuniform cross sec- 
tion, which both introduce vertical and horizontal modes, twist 
produces displacements in al l direc tions. This agrees with the 
results of iRuderman & Scotll (1201 Ih . who considered kink os- 
cillation of a nonplanar coronal loop involving magnetic twist. 
In view of these results the role of twist on stable transverse mo- 
tions has been underestimated in the context of loop oscillations. 
Potentially, the signatures of twist on standing transverse oscil- 
lations could be used as a way to infer the value of the azimuthal 
component of the magnetic field in coronal loops. This could 
have clear seismological applications. However, the fact that real 
coronal loops are in many cases nonplanar and noncircular, i.e ., 
they have an helical shape (see the review o f lAschwandenl2009l) . 
might complicate the comparison between theory and observa- 
tions. 

Finally, it is worth to notice that the analysis performed in 
this paper is based on the assumption of very weak twist. This 
has allowed us to avoid some complications that appear when 
the component of the magnetic field is increased. First of all, 
Alfvenic resonanc es cannot be avoi ded for moderate twist (see 
the work of Karami & Baharill2010h . Second, the modes might 
become kink-unstable. A detailed investigation of these effects 
on standing kink modes need to be addressed in future works. 
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Appendix A: 

The ideal linearised MHD equations in the zero-j6 regime in the 
presence of magnetic twist are the following: 



l(OVr = 



A'oP W oz or J ar 
1 (drb^ . \„ \d{rB^) \ 



lOJVa, - 



lOJV. 



1 
1 



- I —bz -—]B,+ br 

OZ r ar 



iojbr 
i(ob^ 
iojb. 



dz 



dVr im 

B^—+B^ — Vr, 

oz r 



HqP \\ r 



B^ + 



dB, 
dr 



dz ^ dz 



dVr 



dB^ 
'~d7' 



1 / drvr dB. 

- -Bz^ rvr—-^ - BJmv^ + BJmv, 

r\dr dr 



(A.1) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 



where v = (Vr, v^, v^) is the velocity and b = {br, b^, b^) is the 
perturbed magnetic field. If a Fourier analysis is performed in the 
z-direction, the derivatives with respect to z have to be replaced 
by ik. 
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